Two Photon Exchange Contributions to Elastic 
ep Scattering in the Nonlocal Field Formalism 



Pankaj Jain, Satish D. Joglekar and Subhadip Mitra 
Department of Physics, IIT Kanpur - 208016, India 



Abstract 

We construct a nonlocal gauge invariant Lagrangian to model the electromagnetic interaction 
of proton. The Lagrangian includes all allowed operators with dimension up to five. We 
compute the two photon exchange contribution to elastic electron-proton scattering using this 
effective nonlocal Lagrangian. The one loop calculation in this model includes the standard 
box and cross box diagram with the standard on-shell form of the hadron electromagnetic 
vertices. Besides this we find an extra contribution which depends on an unknown constant. 
We use experimentally extracted form factors for our calculation. We find that the correction 
to the reduced cross section is slightly nonlinear as a function of the photon longitudinal 
polarization e. The non-linearity seen is within the experimental error bars of the Rosenbluth 
data. The final result completely explains the difference between the form factor ratio 
Ge/Gm extracted by Rosenbluth separation technique at SLAC and polarization transfer 
technique at JLAB. 
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1 Introduction 



The electromagnetic form factors Fi and F2 parametrize the vertex of electromagnetic in- 
teraction of a photon with an on-shell proton, 

r,b,p') = l,F,{q') + ^/2{q>,.q\ (1) 

where p and p' are the initial and final proton momenta, Mp is the proton mass, Kp its anoma- 
lous magnetic moment and q = p' — p is the momentum transfer. The functions Fi and F2 
are called the Dirac and Pauli form factors respectively. They are experimentally measured 
by elastic scattering of electrons on protons, assuming that the process is dominated by one 
photon exchange diagram (Fig. [1]). We also define = —q^ > 0. Besides the form factors 
Fi and F2, it is also convenient to define the electric and magnetic form factors (or the Sachs 
form factors), Ge and Gm which are more suitable for experimental extraction, 

Ge{Q') = F,{Q') - TKpF^iQ') 

GMiQ') = F,{Q^) + KpF^iQ^) (2) 

where r = Q'^/AM'^. At = 0, Fi = F2 = 1 and Ge = G]\j/fip = 1, where fip is the 
magnetic moment of the proton. The form factor Gm ^ fJ-pGo where is the dipole 
function, 

\^ ^ 0.71/ 

At low momenta, G^; is also approximately equal to Gd- At large momenta, >> 1 GeV^, 

1 

The experimental status of G^; and F2 is, however, currently unclear at large momentum 
transfer. 

A standard technique for the extraction of the proton form factors is the Rosenbluth 
separation [T]. Here one considers the unpolarized elastic scattering of electrons on target 
protons. In the one photon exchange approximation the cross section can be written as 

da auott 



Gm,F,^—. (4) 



dfl e(l + t] 



[tGKQ') + sGliQ')] (5) 



where e = l/[l+2(l+r) tan^(6'e/2)] is the longitudinal polarization of the photon and 6e is the 
electron scattering angle. One finds that the reduced cross section, a^ = TGlj{Q'^)+eG'^{Q'^) 
depends linearly on e. By making a linear fit to the observed a^ as a function e at fixed Q^, 
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Figure 1: The one photon exchange diagram contributing to the elastic electron proton scattering. 
Here k, k' refer to the initial and final electron momenta and p, p' to the initial and final proton 
momenta respectively, q = k — k' = p' — p is the momentum exchanged. 

one can, therefore, extract both Gm and Ge- At large Q^, Gm dominates at all values of 
e. Hence the uncertainty in the extraction of can be large at large Q^. Recent results 
for Rosenbluth separation are available from SLAC [2], [3] and JLAB |1] . The SLAG data 
shows that ^ 1 upto momentum transfer ~ 6 GeV^. The JLAB data is available 



at = 2.64,3.20 and 4.10 GeV^ and shows a similar trend. This result also implies that 
the ratio F2/F1 ocl/Q^. 

A direct extraction of the ratio Ge/Gm is possible by elastic scattering of longitudinally 
polarized electrons on target proton e + p^e + p|5]. In the one-photon exchange approx- 
imation, the recoiling proton acquires only two polarization components. Pi, parallel to the 
proton momentum and Pt, perpendicular to the proton momentum in the scattering plane. 
The ratio. 



where and E'^ are the energies of the initial and final electron. This technique, therefore, 
directly yields the ratio Ge/Gm- The results [HI HIE], available from JLAB, show i^pGe/Gm 
decreases with . A straight line fit to the data gives 



in the momentum range 0.5 < < 5.6 GeV^. The ratio, therefore, becomes as small as 0.2 
at = 5.6 GeV^, the maximum momentum transfer in this experiment. The polarization 
transfer results also imply that QF2/F1 ~ 1 for > 1 GeV^. The observed trend in the po- 
larization transfer experiment is, therefore, completely different from what is measured using 




(6) 




1.06- 0.15 g2 



(7) 
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the Rosenbluth separation. This is clearly a serious problem and has attracted considerable 
attention in the literature [SI [ID] • 

2 Two Photon Exchange 

An obvious source of error is the higher order corrections to the elastic scattering process. A 
reliable extraction of the form factors requires a careful treatment of the radiative corrections 
including the soft photon emission, which give a significant correction to the cross section [Til 
Wli dSl E] . These contributions are calculated by keeping only the leading order terms in the 
soft photon momentum. Furthermore only the infrared divergent terms, which are required to 
cancel the divergences in the soft photon emission, are included in the radiative corrections. 
It is possible that the terms not included in these calculation may be responsible for the 
observed difference. Any such correction is likely to be small and hence cannot significantly 
change the results of the polarization transfer experiment. However a small correction to 
the Rosenbluth separation could imply a large correction to the extracted form factor Ge- 
A possible correction is the two photon exchange diagram which has attracted considerable 
attention in the literature [T^ [TBI \T7\ \TE\ [TU] . Such a diagram is taken into account while 
computing the radiative corrections, but only the infrared divergent contribution is included. 
It is possible that the remaining contribution gives a significant correction. One may also 
consider next to leading order corrections in the soft photon momenta to the soft photon 
emission diagrams. Both of these contributions receive unknown hadronic corrections and 
cannot be calculated in a model independent manner. 

In this paper we estimate the two photon exchange contribution using an effective non- 
local Lagrangian. The box and cross-box diagrams which contribute are shown in Fig. EJj^a) 
and[2](b) respectively. As discussed later, in the non-local formalism we need to evaluate one 
more diagram. The two photon contribution has also been obtained by model calculations 
in Ref . [T71 [TH] . The authors find that they are able to partially reconcile the discrepancy. 
The results of Ref. [171 HE] show that the predicted Rosenbluth plots are no longer linear 
in e. The experimental results obtained from JLAB [1] show very little deviation from 
linearity. The SLAG results [3] can incorporate some non-linearity due to the presence of 
relatively larger errorbars. The present limit on the deviation from linearity is given in Ref. 
[20] . In Ref. [19] the authors argue, using charge conjugation and crossing symmetry, that 
two photon exchange contribution must necessarily be nonlinear in e. If the two photon 
exchange contribution shows large non-linearity as a function of e then it cannot provide an 
an explanation of the observed anomaly. 
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Figure 2: The two photon exchange diagrams contributing to the elastic electron proton scattering: 
(a) box diagram and (b) cross-box diagram. 

3 General electromagnetic vertex of proton 

The elementary electromagnetic vertex of an on-shell proton is given in eq. [H When the 
proton is off-shell, the vertex is expected to be more general. Further, it must satisfy the WT 
identity, following from gauge-invariance, that implies a relation between r^{p,p') and the 
inverse proton propagator, Sp^{p). A local theory of interaction of a proton and a photon 
would have a U{1) gauge-invariance, implied by local transformations and would imply the 
WT identity: 



This identity would be violated if in calculating the two photon exchange diagrams one uses 
the standard on-shell form factors defined in eq. [H and a free proton propagator. Here we 
are interested in formulating the theory in terms of an effective nonlocal action, which will 
allow us to maintain gauge invariance in the presence of form factors in the electromagnetic 
interaction of proton. It is certainly possible to maintain gauge invariance in a local theory 
also but in this case the form factors will arise only after we take into account loop corrections 
in strong interactions. It is not clear how to systematically do calculations in such a case. 
In the present case the form factors are present at the tree level interaction of photon with 




(8) 
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proton. The vertex r^{p,p') satisfies a generalized non-local version of tlie WT identitjE: 

9 {q') q'T,{p,p') = S^\p') - S^\p). (9) 

wliere g (g^) is a function of appearing in the gauge-transformation equations, ultimately 
to be related to a form-factor in the next section. As we shall see in the next section, this 
identity follows from a non-local electromagnetic invariance, and in fact is more appropriate 
for an extended object like a proton. In the local limit, the function glq"^) — > 1, and the 
identity in eq. reduces to the local WT identity. On account of the charge-conjugation 
invariance of the proton-photon interaction, the vertex r^(p, p'), a 4 x 4 matrix, must satisfjj^ 

C-'T,{p,p')C = -Tl{-p',-p) (10) 

where C is the charge- conjug ate matrix, with C7^C~^ = -7J [22] • We now express the 
vertex in its most general form, employing the 16 linearly independent Dirac matrices: 
1, 75, Iti, 7m75, <^fiu and the 4-vectors: Pf" = {p + p'Y , q^" = [p' - pY. 

^i^ip,p') = aP^ + bq^ + c-f^ 

+ dfP^ + efq^ + f(^P^ + g(^q^ 

+ ka^pP''q^P^ + la^pP''q^q^ + m^^-f,C;'^Ppq,. (11) 

Here, the 12 coefficients a,b,....,m are functions of the three Lorentz invariants p'^,p''^,q^. 
Charge conjugation requires that a, c, d, g, j, k are symmetric under p'^ ^ p'"^ and b, e, /, h, I, m 
are antisymmetric under the same operation. To implement the WT identity, we express, 

S^\p) = a {p') p + P {p') (12) 

in its most general form. We then impose the WT identity given in eq. M This leads to 
some constraints between the coefficients. The net result of all this is to yield the following 
form for F^: 

r^(p, p') = a'Pfj, + c'7^ -|- ja^aO^ + d'fP^ + 7 divergence free terms. 



^Such non-local WT identities generally occur in non-local quantum field theories. See e.g. Ref. [21] . 
This WT identity reduces to the usual one as q Q, provided 17(0) = 1. 

^The negative signs for momenta on the right-hand-side are a consequence of our different sign conven- 
tion regarding the incoming particle (incoming momentum positive) and the outgoing particle (out-going 
momentum positive). 
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We enumerate the divergence free (i.e. with g^X^ = 0) terms: 

h' [(p'2 - p')q^ - g2p^] + / (-P.g7^ + ^P^) + g {-q^^^ + ^g^) 
+ k {-a^^P'^P.q + a^pP'^ql'P^) + / (a^«P"g' - a^.P'g^g^) 

We further note the relations that arise from the WT identity and which restrict the form 
of some of the coefficients (a', c') considerably: 



whereas the coefficients of j, 6, /, 5^, /c, /, m are completely arbitrary functions of the Lorentz 
invariants. We make several observations: 

1. We note first that power counting would associate all operators except those three with 
coefficients a', c', j with a local operator of dimension 6 or higher. 

2. We note that the dependence on q^ of both a' and d are identical. Near mass- shelil, 

{[3r + 0[p^-Ml]]g~^ (g^); 

3. The on-shell expression (eq. [T]) for r^(p,p') takes operators of dimensions 4 (electric) 
and 5 (magnetic) into account. It is then logical that the only other operator of 
dimensions 5 should also be included in the off- shell expression for the r^(p,p'). We 
shall take these three terms into account in our minimal effective Lagrangian model. 



a (p2) ^ ao + - M^], [3 (p^) ~ /3o + Pi[p^ - M^]; and thus, 



a' 



4 Effective Lagrangian Model 

We represent the interaction of the photon-proton system by an effective nonlocal Lagrangian 
model based on the discussion in the last section. We adopt the following guidelines in the 
process: 

•^thc condition that Sf{p) ^ f~hr ^^^ar mass-shell requires that q;qMp + /3o = 0; ao + 2Mpai + 2Mpj3i = 1 



7 



The Lagrangian model should incorporate up to dimension 5 operators, for reasons 
partly explained in the previous section. The assumption is that in the effective La- 
grangian approach, the higher dimension operators will contribute much less. This is 
borne out in the calculations performed. (See Figure 15 and the subsequent discussion.) 

The model should incorporate the results regarding the form of the coefficients a', c' 
obtained earlier (see eq. [T3|) . thus at least embody the form-factors on mass-shell. The 
resulting model is necessarily non-local. 

We assume that the model has lowest order derivatives of fermions. Our assumption 
about the dimensionality of operators is consistent with this. 

We require that this non-local model has an equivalent form of gauge-invariance. Such 
constructions of non-local versions of local symmetries are known in literature EI] 
and we shall show explicitly that our model below has a very simple form of non-local 
gauge-invariance . 



A Lagrangian model which satisfies these constraints is given by. 



C = 'ip 



ef[ 



A2 



92 
A2 



b" 
2M 



where iD = id — ef[ ^ A is the non-local covariant derivative. We point out that the 
form factors, f[ and are to be extracted directly from experiments. We make a number 
of observations regarding this effective Lagrangian: 

1. £ is invariant under the non-local form of gauge transformations: 



or equivalently, 



/-I 



A2 



In the latter form, the gauge transformations are similar to the usual local ones, with 
the exception that in the first of these d^f[~^ ff^)- This leads to the non- 

„2 ' 



m 



local WT identity of eq. [9l i.e., with a replacement — > q^j,f]~^ 
eq. [HI Under this transformation, F^^^ and hence the second term is gauge-invariant 
independent of the form of Also, the non-local gauge-covariant derivative satisfies: 
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2. The last term generates a term proportional to in r^(p, p') with a form factor 
proportional to f[, the same one appearing in the electric term. This is consistent 
with the comment on the form of a' and c' given earlier. 

3. The (non-local) gauge- invariance of the last term requires that it is composed of the 
(non-local) gauge-covariant derivative: this restricts the form-factor present in this 
term as above. 

4. The Lagrangian is exactly valid as long as the proton is on-shell, irrespective of the 
value of the momentum transfer g^. In this limit the interaction of the proton with 
photon is described in terms of the two form factors and no higher dimensional terms 
are required. The higher derivative terms we drop give higher order contributions in 
powers of (P^ — ) / A^, where P^ — is the offshellness of proton momentum. These 
higher order terms can be dropped as long as the dominant contribution to a process 
is obtained from the kinematic region where (P^ — M^)/ « 1. 

5. In the limit A — ^ oo we reproduce the local field theory model for a proton with an 
anomalous magnetic moment. 

Before we proceed, a comment on the non-local form of gauge-invariance is in order. It 
appears that a local form of gauge-invariance for extended particles such as a proton is 
inappropriate. Consider the wave-function of an extended particle centered at x: viz. tp{'x,t). 
Let y be a point within the charge-radius R of the proton: |x — y| < R. Let us imagine that 
a gauge-transformation on is carried out (at t) around y with a very narrow support, 
p : p << |x — y|. In the model of fundamental constituents, the quark wave-function should 
be affected around y, which in turn should affect the proton wave-function even though the 
gauge-transformation at x, depending on a (x) will be zero. Thus, the proton wave-function 
should be affected by a local gauge-transformation with a support anywhere in its charge 
radius. The above form of non-local version of gauge-transformations embodies this idea. 
Note that the Fourier transform of f[ ^ has a support over a distance ~ 1/A ~ P. 

It proves convenient to rearrange the Lagrangian as followj^ (Recall the relation = 



b - /.cz _\2 



^Actually, the constant Alp and the normahzation of KE term are also modified below. However, we shall 
soon modify the form of the Lagrangian further, where this proves unnecessary. 
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Had there been no magnetic term, the last term would have formally vanished by classical 
equation of motion. We note that an inclusion of the last term has now modified the inverse 
propagator: it has non- vanishing terms at e = 0. This, in particular, gives a spurious pole 
in the propagator at another value of p. This problem can be avoided if we can write the C 
in the following form: 

M,) e.p {,3 - M,) } ^ + [^] i, (14) 

which is now understood to have been consistently truncated to a given order in b. We now 
note that the inverse propagator is: 



and has only one zero at p—Mp = and the residue of the propagator at the pole is 1. (In this 
form of C, a is related to the anomalous magnetic moment, Hp by the relation, d = eKp/2 and 
Mp is the physical mass). Since we shall consider, in the two photon exchange calculation, 
terms with the last operator of dimension 5 inserted in two photon vertices, we shall do the 
entire calculation consistently to O (6^) using eq. [TH In this case, the propagator for the 
proton is 

FiJ^ {'A ^* ~ " - ^ + 

5 Reduction of the action 

In this section, we shall find an effective way to calculate the matrix elements involving 
insertion of the last term in the action. Since a 2-photon exchange diagram at 1-loop is at 
most we shall evaluate the effect of this term to . What we are interested in are 

the tree-order matrix elements of two (possibly virtual) photon emission from an on-shell 
proton. The calculation of these can be simplified considerably in this context with the use 
of the fermion equations of motion. The result is simple: of all the terms up to viz. 

6a, P, 6^a, Pa^], only the last term of (9[Pa^] gives a non-zero result. While the result 
can be worked out, it is most effectively dealt with in the path-integral formulation. 
We define, 

W[J^, K,K] = J D4>exp i I j d^x[C + J^A^ + Kij+ ^/sfj 
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where L is the action of eq. [TH and denotes generically the measure of the path-integral. 
We now perform a field transformation: 



■0 = exp 



2M„ 



i1) - Mp) 1 



Under this transformation, the Jacobian is, 



J = det exp 



2M. 



(^i3 - My 



exp tr 



2AL 



W - Mr 



exptr -X ^ = a constant. 



and hence can be ignored for the connected Green's functions. This then yields. 



expi< / (i^x 



X exp 



+ i^exp 



2M„ 



2Mp^ V 
X expi<v / (i^x 



92 



A2 



2M„ 



Mr, 



(15) 



(16) 



We note that if we had a = 0, we would have no left-over term in the action. It is easy 



to show that K 



exp 



2Mr, 



proton matrix elemental- Thus all tree- 



■ip' does not contribute to tree-level on-shell 
evel on-shell 2-proton matrix elements involving the 



last term in eq. [THare at least of order ab. We now expand the action to We find, 

WiJ", K,K] = j D(P expiS^j d^x [C + J^A^ + K^j' + i^K] | + 0[P] 



with 



+ 2^^ (^-^^ 



+ 



2Mp V 









_A2_ 









_A2_ 




1^' 



2Mp 8M2 



^This is because one does not have a pole in at least one of the external momenta: it is forbidden either 
by an explicit factor of ^ — Mp, or by a vertex. 
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We now perform another field transformation, 



A2 



+ 



2Mp 8M2 



it) - Mr 



+ (17) 



We can write, 



V' = ^ i 1 + 7Tirr^M^/2 



2M„ 









_A2_ 



2Mp 8M2 



Mr 



2M. 



d 



2 H 



A2 



2Mp 8M2 



The (9[F2] will not matter for the present calculation of 2-photon exchange, as it will give a 
term having 3-photon fields. The Jacobean for this transformation. 



I /J' = det[l + X] = det 
" ■ X2 



Y2 v2 



det 



exp [trX] det 



det det 



1 - e 



-X 



1 - e 



X2 



1 — tr I e 



-X 



+ C[P 



X2 

1 - tr— + 0[b^] 



32M4 



A2 



A2 



F^'' X (constant) 



The last term does not contribute to the emission of two photons from a proton line in the 
tree approximation. As a result of the transformation (eq. [T71) . the action then becomes. 



2M, 



2M„ 



2M 



fe2 \ 



A2 
A2 



A2 



UjD-M,) la, J!, 



A2 



and the source term transforms into 

a 



4j'[i + x]-^k 



^ 1 + 



— / / a 

+ ip 



2Mr, 



2Mp 



'^fj,iyf2 



A2 



2Mp ^ 8M2 



Mr 



K 



A2 



A2 
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Figure 3: The two photon exchange diagram proportional to b'^ contributing to the elastic electron 
proton scattering. The Feynman rule for this diagram can be obtained from eq. [TSl 



None of these terms contribute to the tree approximation 2-photon matrix element for reasons 
similar as before. 

In conclusion, when we look at the 2-photon exchange diagrams having up to two in- 
sertions of the last term in eq. [TU each set of diagrams contain a common part, viz., two 
(unphysical)-photon tree amplitude from an on-shell proton. The above discussion shows 
that the net effect of that comes from the terms 



2M„ 



2M„ 



b \-, 



2M„ 



Mr 



ij' (18) 



We shall show in appendix 1 that the first term does not contribute in the Feynman gauge. 
That leaves us with only the last term. The Feynman diagram corresponding to this term 
is shown in Fig. [3l 



6 Calculation and Results 

In this section we give details of the calculation of the two photon exchange diagrams using 
our effective Lagrangian. The calculation turns out to be complicated due to the explicit 
presence of form factors at the vertices. We also require models for the form factors both in 
the space-like and time-like regions. In the space like region the form factor Fi{q^) is known 
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Figure 4: The amplitude of Gm/ jJ^p (solid line) and Ge (dotted line - Model I, dashed line - Model 
II). The unfilled squares represent the constrained values for G^; at = and = 4M^. Results 
of Rosenbluth extraction experiments (filled circles - JLAB, unfilled circles - SLAC) are also shown. 



reasonably well. In the time like region experimental data exists for the form factor Gj\/(g^) 
for 4Mp < < 14 GeV^, where 4Mp is the threshold energy for pp production. In Ref. 
[21], Gm{(i'^) has been extracted in the unphysical region < < 4Mp by using dispersion 
relations The extracted form factor shows two resonances at masses M ~ 770 MeV 

and M ~ 1600 MeV. The phase of the magnetic form factor also shows a large variation in 
the unphysical region. The electric form factor GE{<f)i however, is not well known. The 
amplitude in the unphysical region is obtained in Ref. [27] • However the phase is not known. 
Our model for the form factors is given in Appendix 2. We use two different models. Both 
consist of a sum of simple poles. The corresponding masses and widths are given in Tables [3] 
and m The values of these parameters are obtained by fits to the experimental data, or the 
data obtained from experiments by using dispersion relations [23]. The resulting amplitude 
and phase of the form factors for the two models are shown in Fig. H] and respectively. 

Using these models for the magnetic and electric form factors we can obtain the form 
factors Fi and -F2, required for our calculation. The models used are convenient since they 
allow us to use the Feynman parametrization to compute the loop integrals. The form factors 
for the two models have a small imaginary part even for space like momenta. However this 
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Figure 5: The phase (in degrees) of Gm (solid line) and Ge (dotted line - Model I, dashed line - 
Model II). 



region contributes negligibly to the loop integrals. The dominant contribution comes from 
the unphysical region < < 4Mp where the form factor is several orders of magnitude 
larger than its value in the space like region. In this region our model provides a very good 
fit to the extracted form factor [21] . Moreover the imaginary part in space like region is very 
small and unlikely to affect our results significantly. The resulting amplitudes and phases 
for model I are shown in Figs. E] and [7] respectively. The corresponding results for model II 
are shown in Figs. [H]and[ni 

Using the form factor models in Appendix 2 we can determine the amplitudes of the box 
and cross box diagram as well as the amplitude due to the extra term displayed in eq. [13 
The box diagram amplitude can be written as, 



iM 



B 



E 

a,b 



(27r)4 {{k - If -ml + - /i2 + i^){q^ - + i^) 



X 



2M„ 



(p + /)2 - M2 + 



2M 



(19) 
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16 





17 



where q = q — I, nie is the mass of the electron and ^ is an infinitesimal positive parameteiij. 
A small mass of the photon fi has been introduced in order to regulate the infrared divergence 
in these integrals. The infrared divergent part has to be subtracted from our result since it 
is included in the standard radiative corrections which are applied while extracting the form 
factor. Using the form factor model in Appendix 2, we find, 



iM} 



E 

a,b 



dH 



X {{C,g,{l) - mlD,g,{l))^, + i2MpD,g,{l)a,^lP}U{p)\ . 



(20) 



It is convenient to rewrite this expression in terms of the coefficients C and D', defined in 
Appendix 2. We find, 



iMi 



E 



dH 



u{k')Y{^- l)Yu{k) 



[U{p'){{C[gS) - ^M^D'Mq)h, + i2M,D'MqKa{qr} 

{{C'^gjil) - AMlD'^g,{l))^, + t2M,D'^g,{l)a,pf}U{p)\ 
dH 1 



{p + ly - M2 + 



= e 



(27r)4 {{k - ly -ml + tOiip + ly - M2 + i^) 
[ni{Cl A', D'^) + A', C'^. D'^)] 



(21) 



(22) 



h3 



where the last step defines the factors ni{C[, D[, Cj, D'-) and n2{C[, D[, Cj, D'-). We may can- 
cel the /2 factor multiplying n2{C-, D'^, Cj, Dj) with a factor {p — + in the denominator. 
We then find, 



iM 



B 



dH 



{2ny {{k - ly -mi + ii){{p + ly - + zo 

[ni{C[D[, C'^. D'^) + yn,{Cl D'^) 

id 



(23) 



^Here we use the notation ^ instead of the standard notation e to avoid confusion with the symbol e used 
to denote the photon longitudinal polarization 
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The cross-box diagram amplitude can be written as, 



iM 



CB 



E 



(2vr) 



{k — ly — ml + 



1 



X 



+ ^ - / + Afp 



Kir. 



a,b 



(2^ 



M(A:')7''(^-/)7^M(fc) 



1 



(A; — /)2 — ml + 

f + (^-t + Mp 
{p + g)2 - M2 + 

{{Ca9a{q) - AM^Da9a{q)h>. + t2MpDaga{qMqr}U{p)] ■ 



(24) 



The amphtude proportional to h"^ is given by: 



8M2 

^ ' a 



E 



(2^ 



{k — ly — ml + iC, 



{P-fi^ + tO{q^-f^' + ^0. 



X U{p') 



IK 



2M, 



%K, 



p 



2Mp 
8M2 



F2{l)a.,l^ ) + ^ - / - Afp) ( ^F2(g)a,„(g)" 



p) 



liKjr 



2M, 



U{p) 



(2^ 



u{k')-i^'{]/i - t)Yu{k) 



X 



(A: - /)2 - m2 + 2^ 

U{p') [{i2MpD[gi{q)a^^{qr){p + / - Mp){t2MpD'^g,{l)a,pl'') 
{t2MpD'^gj{l)a,fsl'')if + ^ - / - Mp)(^2MpD:(7,(g)a^,(g)")] 



(25) 

In our numerical calculation we set the mass of the electron m^ = 0. 

The contribution of the two photon exchange diagrams to the electron-proton elastic 
scattering cross section can be written as 



dVtp 



+ 0(a'). 



where, 



= Mb^ Mcb + M 



6) 



(26) 



(27) 
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is the total amplitude of the two photon exchange diagrams and 

Mo = -^u{k')Yu{k)U{p') (F,{qh, + U{p) (28) 

the tree amplitude. Hence, the contribution of the two photon exchange diagrams to the 
reduced cross section is given by: 
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e{l + t) da^^ 

O'Mott dVte 



^4, 



g'£(l + T) 



The diagram proportional has no infrared (IR) divergent term. So the contribution 
coming from it is computed keeping /i^ = 0. Contribution from box and cross-box diagrams 
are computed at 10 different values of /j,^ (from 0.005 to 0.0095). The numerical calculation 
of the cross-box diagram is straightforward since the integral is well defined. However for 
the evaluation of the box diagram the numerical evaluation is facilitated by keeping a small 
imaginary term ^ in the propagators. This makes the integral in the infrared limit well 
defined in the case mg = 0. For each value of g^, e and /i^ we have calculated the box diagram 
amplitude for 4 different values of ^ (between 0.001 and 0.00175). The amplitudes depends 
almost linearly with ^. The final dependent box diagram amplitudes are obtained by 
extrapolation to ^ = 0. The two different models for the form factors described in Appendix 
2 gives almost identical results. So for the rest of the section we quote the result obtained 
using Model-I only. 

The IR behaviour of the two photon diagrams have been calculated by Mo and Tsai 
p!2l [T3] . In the limit /x^ — the leading term from the box and cross-box diagram can be 
expressed as: 



M'jI = -[K{p\ k) - K{p, k)]Mo, (30) 



a 



TX 



where 



dx 



[xpi + (1 - x)pj)'^ 



The IR contribution to the reduced cross section coming from the box and cross-box diagram 
is given by: 

a'jl=^-^[K{p',k)-K{p,k)]ai. (31) 
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Let 



cr 



.27 = 
IR — 



dir + bir In jj}. 



(32) 



To remove the IR part from we fit it with the following function: 



.27 

R 




(33) 



with 



Here gives the IR removed cr^'^. It has been explicitly verified that keeping 0[fi^] terms 
in a(/i) and b{fi) has no effect on the slope (thus on Ge) of cr'^{e) with respect to e. The 
difference between these two fits leads to a very small correction to Gm only and hence can 
be ignored. 

The result of the calculation for the box and cross-box diagrams is given in Figs. [TUland 
[TTl Here we have considered the momentum transfer = 2.64, 3.20, 4.10, 5.00, 6.00 GeV^. 
The first three values are same as those used in the JLAB extraction of form factors using 
Rosenbluth separation. The contribution from the diagram proportional to 1? is shown in 
Figs. [12] and [131 Here b is taken as 1. We fit cr^^^ = (cr^ + cr'^^) and cr|j (here and for rest 
of the section we use the notations a'^ , to denote the IR removed contributions) to 

the following functions: 



From Table [T] and [2] we also see that the contribution due to the b term is relatively small 
as long as the magnitude of b is of order unity. As the magnitude of b is unknown we shall 
assume ~ and take ~ ^r the rest of the section. 

Fig. [TH shows the contribution of the dimension five operator proportional to F2 to the 
reduced cross section. This contribution is obtained from the box and cross-box diagrams. 
For comparison we also show the total contribution of both these diagrams. The IR ^ 
dependence is not removed in this calculation and the parameters chosen are = 4.10 
GeV^, = 0.005 GeV^. We find that the contribution from terms proportional to F2 x F2 
is much smaller compared to the total contribution, justifying the truncation of our action 
to only operators of dimension 5. 



hie) 
f2{e) 



ci + C2e 

di + d2e + d-iS^ 



(34) 
(35) 



The values of Ci, C2 and c/i, c^s for o"^*"^ and are given in Table [T] and [2] respectively. 
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Cl 


C2 


di 


d2 


^3 


2.64 


-3.39 
±0.30 


5.59 
±0.48 


-4.06 
±0.12 


9.58 
±0.59 


-4.04 
±0.58 


3.20 


-2.31 
±0.25 


3.99 
±0.45 


-2.94 
±0.08 


7.53 
±0.38 


-3.69 
±0.38 


4.10 


-1.20 
±0.15 


2.23 
±0.25 


-1.61 
±0.11 


4.30 
±0.47 


-1.95 
±0.43 


5.00 


-0.79 
±0.16 


1.63 
±0.28 


-1.13 
±0.10 


3.59 
±0.47 


-1.96 
±0.46 


6.00 


-0.44 
±0.12 


1.10 
±0.21 


-0.69 
±0.08 


2.55 
±0.40 


-1.45 
±0.39 



Table 1: Values of fitting parameters for (cj^ + ^r^)- numbers have been scaled by 10^. The 
parameters ci and C2 are defined at eqn. [3lland di, d2 and da are defined at eqn. [35l The standard 
errors obtained in fitting the result are also shown. Figs. (TU] and [TT] shows the functions /i and /2 
for this case. 





Cl 


C2 


di 


d2 


dz 


2.64 


11.60 
±0.86 


-8.23 
±1.41 


9.65 
±0.37 


3.37 
±1.81 


-11.73 
±1.79 


3.20 


10.21 
±0.45 


-5.27 
±0.83 


9.04 
±0.12 


1.29 
±0.59 


-6.84 
±0.60 


4.10 


8.75 
±0.31 


-3.67 
±0.52 


7.86 
±0.04 


0.86 
±0.16 


-4.27 
±0.14 


5.00 


7.37 
±0.19 


-2.23 
±0.34 


6.95 

±0.06 


0.20 
±0.29 


-2.43 
±0.28 


6.00 


6.36 
±0.12 


-1.55 
±0.21 


6.12 
±0.09 


-0.12 
±0.44 


-1.43 
±0.42 



Table 2: Values of fitting parameters for o"^ with 6 = 1. All numbers have been scaled by 10^. 
The standard errors obtained in fitting the result are also shown. The parameters ci and C2 are 
defined at eqn. [HUand di, d2 and are defined at eqn. ESI Figs. WI\ and [TH] shows the functions 
/i and /2 for this case. 



22 



(a) q2 = 2.64 GeV^ 
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(b) = 3.20 GeV^ 




(c) Q^ = 4.10Gev2 
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Figure 10: Total contribution of the box and cross-box diagram to the elastic electron proton 
scattering for = 2.64 GeV^ (a), 3.20 GeV^ (b), 4.10 GeV^ (c). The dashed lines represent /i 
(eqn. and the solid curves represent /2 (eqn. [35|) . The fitting parameters are given in Tabled) 
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(b) = 6.00 GeV^ 
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Figure 11: Total contribution of the box and cross-box diagram to the elastic electron proton 
scattering for = 5 GeV^ (a), 6 GeV^ (b). The dashed lines represent /i (eqn. and the solid 
curves represent /2 (eqn. [35]) . The fitting parameters are given in Table [TJ 
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Figure 12: Contribution of the diagram proportional to 6^ to the elastic electron proton scattering 
for Q2 = 2.64 GeV2 (a), 3.20 GeV^ (b), 4.10 GeV^ (c). Here b = 1. The dashed lines represent /i 
(eqn. and the solid curves represent /2 (eqn. [35|) . The fitting parameters are given in Table [2j 
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Figure 13: Contribution of the diagram proportional to W' to the elastic electron proton scattering 
for = 5 GeV^ (a), 6 GeV^ (b). Here 6=1. The dashed hues represent /i (eqn. [M]) and the 
solid curves represent /2 (eqn. [35]) . The fitting parameters are given in Table [2j 
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e 

Figure 14: The contribution to the reduced cross section coming from terms proportional to 
F2{l)F2{q - 1) in McB and Ms for = 4.10 GeV^, //^ = o.005 GeV^ (filled circles). The unfilled 
circles represent the total contribution coming from both the form factors Fi and F2 for the same 
Q'^ and /i^. 

To obtain the corrected cr/? we subtract the linear fit to a^^ , fVi.^) (see eqn. [311 Table 
[1]), from a linear fit to a^{e) given by: 

(7^R{e) = Q^{l + eGi). (36) 

Then the corrected reduced cross section, 

cjnie) ^ 4^ - = {g, - ) + 5(^0^1 - c^). (37) 

In Fig. dSlwe plot au^e) for different . We determine the corrected form factors Gm 
and Ge by: 

Gm = \\IQo-c;' (38) 

Ge = \Iq^Qx-c^ (39) 

Fig. [m shows how the ratio ix^GeIGm is modified by the two photon exchange contributions. 
The SLAG Rosenbluth data after applying the two photon exchange correction is shown by 
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(b) q2 = 3.20 GeV^ 




(c) Q^ = 4.10Gev2 




Figure 15: Corrected cross section, aji (solid lines) obtained using eqn. EZlwith 6 = for Q2 = 2.64 
GeV^ (a), 3.20 GeV^ (b), 4.10 GeV^ (c). The unfilled circles represent the data points obtained by 
Rosenbluth separation method at JLAB and the dashed lines are the straight line fits to these. 
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Figure 16: The ratio, fJ-pGE/GM obtained by polarization transfer technique at JLAB (filled 
triangles) and Rosenbluth separation technique at SLAC (unfilled squares) and JLAB (filled squres). 
The ratio after correcting for the two photon exchange contribution is also shown. The filled circles 
are the corrected JLAB Rosenbluth data and the dotted line is the best fit through these points. 
The unfilled circles are the corrected SLAC Rosenbluth data and the dashed is line the best fit 
through these points. 



the unfilled circles. The dotted line represents the best linear fit through this data. We 
find that the two photon exchange correction completely explains the difference between the 
SLAC Rosenbluth separation data and the JLAB polarization transfer data. However it is not 
able to explain the difference between the JLAB Rosenbluth and polarization transfer results. 
The corrected JLAB Rosenbluth data is shown by filled circles. The JLAB Rosenbluth data 
lies systematically above the SLAC data. 



7 Conclusions 

In this paper we have constructed a nonlocal Lagrangian to model the electromagnetic inter- 
action of proton. The model is invariant under a nonlocal form of gauge transformations and 
incorporates all operators up to dimension five. The model displays the standard electro- 
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magnetic vertex of an on-shell proton. The dimension five operators also contain an operator 
with an unknown coefficient whose value can be extracted experimentally. We use this model 
to compute the two photon exchange diagrams contributing to elastic scattering of electron 
with proton. The calculation requires the proton form factors in the entire kinematic range. 
We find that the two photon exchange diagram contribution to the reduced cross section 
aR shows a slightly non-linear dependence on the longitudinal polarization of the photon e. 
The non-linearity seen is within the experimental error bars of the Rosenbluth data. We 
apply the correction due to two photon exchange contributions to both the SLAC and JLAB 
Rosenbluth separation data. The resulting cross section for the SLAC data is completely 
consistent with the JLAB polarization transfer results. However the JLAB Rosenbluth data 
still shows a large deviation. It, therefore, appears that the two photon exchange is able to 
explain the difference in the experimental extraction of proton electromagnetic form factor 
Ge using the Rosenbluth separation and polarization transfer techniques if we accept the 
SLAC Rosenbluth data, which is available over a larger momentum range. 

8 Appendices 
8.1 Appendix 1 

We shall show in this appendix that the first term in eq. [18] proportional to 

does not contribute to the 2-photon matrix element in the 1-loop approximation in the 
Feynman gauge in the zero electron mass limit rrif. = 0. For this purpose we write the term 
as 

We then note that 

• o-^^cTAp + o-ApO-^.;. = a linear combination of I and 75 = a {g^xOup - 9p,p9u\) I + ip^p,^Xp75, 
where a, (3 are constants and, in particular, no a— terms appear. 

• The Feynman integral has no dependence on both p and p' . 

• Thus, the result for the 2-photon exchange diagram is of the form: 

u{k')'y^'yaluuik)U{p')[l^,]U{p) X I^^^^{k,k') 



A2 



2 



A2 
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On simplification, this becomes 

+ 42eMa./37^5} u{k)U{p') [I, 75]t^ {p) X I^^^{k, k') 

Now, g^al'^°"^ik, k') is a linear combination of terms that are ~ k'^ or k"^. Both of these 
terms give zero. Similar logic applies to all other terms. 



8.2 Appendix 2: Model for the Form Factors 

The fits for Ga///^p and Ge are given by: 

.2^ 4 



^ (g2 -ml + imaV'^) 
B' 



(40) 
(41) 



We have considered two fits for Ge- The values of the masses and the parameters are 
tabulated in Table (Model I) and Table 1 (Model II). 



a 


K 


b: 






1 


-2.882564 
+z 1.944314 


-3.177877 
+z 2.123389 


0.8084 


0.2226 


2 


2.882564 
-z 1.944314 


3.177877 
-z 2.123389 


0.9116 


0.1974 


3 


-1.064011 
-z 3.216318 


-0.608148 

-z 5.685885 


1.274 


0.5712 


4 


1.064011 

+z 3.216318 


0.608148 
+z 5.685885 


1.326 


0.5488 


5 





3.211388 
+z 0.693412 


1.96 


1.02 


6 





-z 0.693412 
-z 0.693412 


2.04 


0.98 



Table 3: Masses, widths and parameter values for Gm/i^-p and Ge fits (Model I). A's and B s are 
defined in eq. HO] and HH 
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a 




b: 


ma 




1 


-2.882564 
944314 


-3.392256 
+z 2.194129 


0.8084 


0.2226 


2 


2.882564 
-z 1.944314 


3.392256 
-z 2.194129 


0.9116 


0.1974 


3 


-1.064011 
-z 3.216318 


1.224037 
-z 6.877523 


1.274 


0.5712 


4 


1.064011 
+z 3.216318 


-1.224037 
+z 6.877523 


1.326 


0.5488 


5 





1.645805 
+z 1.824298 


2.107 


0.663 


6 





-1.645805 
1.824298 


2.193 


0.637 



Table 4: Masses, widths and parameter values for Gm/ l^p and Ge fits (Model II). A s and B s are 
defined in eq. SO] and HH 

Using the models for the magnetic and electric form factors we can determine the Dirac 
and Pauli form factors. Let the fits to the form factors Gm and Ge be: 



a=l 
6 

GE{q') = Y.Baga{q' 



(42) 
(43) 



a=l 



with Aa = fJ-pA'^ and Ba = B'^. The Qa^ are defined by, 

9a{q 



(44) 



The form factors Fi and F2 are given by, 

Gm — Ge . , ,2 — Gm 
HpF2 = — : = 4M, 



1 + T '-'P q^- 4M2 

6 

Fi = Gm - KpF2 = Y^iGaQa - ^M^Da~9a) 



(45) 
(46) 



a=l 
-.2 A T\l2\ 



where Ga = Aa, Da = Ba — Aa and Qa = Qa/ (q — 4Mp). These definitions are convenient in 
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evaluating the two photon exchange amphtudes. We also have, 

6 



E 



g2 _ ^2 + ^ (g' - + ^0 - 4M2) (g2 - m2 + zFj 



E 



with my = /i, mg = 2Mp, Fy = ^ and Fg = 0. Here (a = 1, 2..., 6) 

1 1 



(47) 
(48) 



D' 



(4M2 - /i2 + 2^) [m2 - 4M2 - tVa ml-fi^ + t^- tTa_ 



= E 



^ (4M2 - /i2 + ^e) (ml -li^ + ii- iFj ' 



E (4M2 - /i2 + ^e) {ml - 4M2 - zF,) 



0. 



The coefficient D'^ is found to be zero and hence the summation in eq. HE] terminates at 
i = 7. Similarly, 

^ C' - AM^D' 



E 



g2 - ^2 _^ ^ g2 _ ^2 _^ . ' 



(49) 



where 



a 



Ca 



ml- iJ? + - iVa ' 

6 



-E 



ml - + iC, - iTa ' 



a = 0. 



We can also write Fi and k„-F2 using this general notation. We find 



K.pF2 



Fi 



E 



4M2 D" 



^ g2 _ ^2 _^ . 
^ g2 _ ™2 _^ 



(50) 
(51) 
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with 

D. 



D'^ = 0, 
C' = A„ 



'■ai 



a; = 0. 

8.3 Appendix 3: Sample Calculation: Box Diagram 

Here we present a sample calculation of one of the terms in the Box diagram. The contribu- 
tion of the box diagram amplitude to the two photon exchange cross section is proportional 
to, 

^^^=v^^ ^^^^ 

id 



where, 



1 



(/2 - m] + iV,) 

We can now evaluate this integral by the standard Feynman parametrization technique. We 
define 

V = f + 2l.{x2P - xik - xsq) + X3(g^ - m^) - X4m^ + + iix^T^ + XiT'j) (53) 

with F- = Fj — C,- We now define the shifted momentum, r = I + {x2P — xik — x^q) which 
gives, T> = r"^ — A', with 

\ f '2, '2 '2, '2, / \2 2 2 

A = X;^mg + XgMp — X3(l — Xi — X2 — Xsjg — 2xiX2EMp + x^m^ + Xirrij 
- i^- i{x3T[ + x^T'j) 

With this momentum shift the numerator becomes: 

=^^o + r^W^ + r^rM^" + r,,r,r^r- 

Hence, 



Xg = 6 j\UdxjC£ 



X. 



a 



(27r)4 (r2 - A 
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As the denominator depends only on the magnitude of r, 



(27r)4 p4 0' 
d^r r^rMr f d'r {g^M^'r^ 



(27r)4 p4 J (27r)4 



d^r r^r^TpM^"^ 



0. 



(27r)4 p4 
Let Af2 be the shorthand notation for ^g^yAfl^^ . 

^ " d^r ATp + rWa 

(27r)4 (r2 - A')4 



With Jo = / dx3 / dx, I dx^^^ 

A 



^0 ^0 
and T2 = / (ix3 / dx2 / rfxi — . 



^0 



Here 



A = X;^mg + XaMp — X3(l — xi — X2 — xsjg — 2xiX2-E'Mp + xsm^ 
+ (1 - xi - X2 - a;3)m^ - - «(a;3r- + (1 - Xi - X2 - Xz)V'-). 

If we neglect the mass of electron then, 

A ^ X2Mp — X3(l — xi — X2 — xsjg — 2xiX2-EMp + X3mj 

+ (1 - xi - X2 - X3)m| -ii- «(x3r^ + (1 - Xi - X2 - x^)^'^) 
= Xxi + Y, 

where, 

X = -2x2EMp + x^q^ - + iV'^ 

Y = xjMp - xs{l - X2 - xs)q'^ + X3{m'^ - m^j) + {1 - X2)m'^ 

- ^e-^x3(^:-^;.)-^(l-x2)^;.. 

A/i, A/2 can be written as: 

A/q = Z-i + Z4X1 + Z^xj 

J\f2 = Z1 + Z2X1. 
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Then 



Jo = dx-i I dx2 I dxi — — — — — 

io Jo [Xxi + Yf 

and ±2 = f dx-i / "3;2 / dxi - 



>/o 



'{Xx^ + YY 



The Xi integration can be done analytically to obtain: 



= L ^"^3 / dX2 



1 r^-x3 r_2ZiX2 + X(2Z2F + Z4) -2Z5F , fXL + Y 

m ' 



Jo 



+ Y) 



X3 V ^ 

{2^5^' + X\Zi - 2Z2LY) - XY{2Z2Y + Z4 - Z5L)} 



(54) 



where L = 1—X3~X2- Zi^s are obtained using FORM [28] and J]^'s are numerically computed 
using Gauss-Legendre integration technique [29] . 
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